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Abstract. Recent experimental discovery of extended self-similarity (ESS) was one of the 
most interesting developments, enabling precise determination of the scaling exponents of 
fully developed turbulence. Here we show that the ESS is consistent with the Navier-Stokes 
equations, provided the pressure -gradient contributions are expressed in terms of velocity 
differences in the mean field approximation ( Yakhot, Phys.Rev. E63, 026307, (2001)). A 
sufficient condition for extended self-similarity in a general dynamical system is derived. 
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Scaling relations for velocity structure functions in isotropic and homogeneous turbulence 
are defined as : 

S n ,m =< (M(x + r,t)- M (x,t)) n (t;(x + r,t)-t;(x,t)) m >= c nim (£r)^(^)«--^0 n>m (-^, -) 

L f L .f V 

(1) 

where u and v are components of velocity field parallel and perpendicular to the displacement 
vector r, respectively The universality assumption implies that the coefficients c„ im = 0(1), 
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independent of the Reynolds number (dissipation scale n pa LfRe~*). The dissipation rate 
£ = (diUj) 2 = const = 0(1) is equal to the power of external kinetic energy pumping. The 
shape of the structure functions (1) is an assumption , not following any rigorous theory. 
In the inertial range (jjj — > , '- — > oo) the scaling functions 4> n m(j') — > Q>n,m = const, 
independent of the displacement r. 

Both physical and numerical experiments show that the functions n m start deviating 
from the constant inertial range values at r/n pa 10. Since one does not have theoretical 
expressions for (f> n>m , accurate measurements of exponents £ n rrt in a fully developed turbulent 
flow requires an extremely wide range of variation of the displacement r which is possible 
only if the Reynolds number of a flow is huge. This problem is even more severe for numerical 
simulations of turbulence, where usually the wide inertial range is difficult to generate. 
It has been shown in a remarkable paper by Benzi et al [1] that even in the medium (quite 
low, actually) Reynolds number flows, where (1) is hard to observe , the following relation 
(ESS) holds: 

S n ,o(r) = C n0 S mfi (r)^ (2) 

where 8(nm) = I 2 -. It is clear from (1) that if c nm are reynolds number independent, 
then the coefficients C n ^ m in (2) do not depend on the dissipation scale n (Reynolds num- 
ber). Since the range of validity of expression (2) is much wider than that of (1), accurate 
determination of exponents 3(nm) enables one to evaluate the exponents C, nm even in the 
not-too-high Reynolds number flows. Comparison of the exponents calculated this way with 
those measured in extremely high Re flows (j3(nm) ~ (,n,o/£,m,o) was usually extremely good 
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[3]. Since its discovery the relation (2) evolved into a major tool for experimental and nu- 
merical determination of the exponents £ nm [l]-[5]. The definition (2) was introduced and 
tested in Ref. [2]. Since S^o oc r in the inertial range, it is is typically used in application of 
the ESS (2) for analysis of experimental data. It is shown below that extended self-similarity 
(2) can be derived self-consistently from the Navier-Stokes equations. 

It was shown that in a statistically isotropic, homogeneous and incompressible flow governed 
by the Navier-Stokes equations the following equation can be rigorously derived in the limit 
r/Lj — > where the forcing function can be neglected [6] : 



dS 2n ,o . d-1 (d-l)(2ra-l) , , , , , , 

dr + -^—S 2n ,0 S 2n -2,2 = ~(2n - l)V X)2n -2 + (2n - l)vD„,2n-2 (3) 

where 



V X ,2n-2 = (PA*') -^(X))(A M ) 2 ™- 2 (4) 



and 



D u , n = (y 2 u(x>) - W 2 u(x))(^uY (5) 

These relations are exact even in the low- Reynolds- number statistically isotropic and ho- 
mogeneous flows in the range r/Lf — > 0. It is important that -D n)2n (V) = 0(1) and thus, 
vD Ui 2n ~ * as v — > in the inertial range. On the other hand, due to the dissipation 
anomaly, vD u ^ n +\ is finite in the inertial range. To prove the former statement, we consider: 



(2n - l)uD U:2n _ 2 = -(2n - l)(2n - 2)(S U (2) + S u (l))(Au)^ + vd 2 r S 2n - lfi (r) (6) 



where S u = v(du) 2 . The second term in (6) disappears in the inertial range in the limit 
v — > 0. To estimate the first contribution, we write neglecting the subscript u: 

(S(2) + S(l))(Au) 2 ^ < ^(S(1)+S(2)) 2 S^ (r) (7) 
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Since in the inertial range (r — > 0), (£(1) + £(2)) 2 oc with /i 0.2, this term is negligibly 
small compared to the O(S 2nj0 /r) contributions to (3) for not too small moment number n, 
provided C, 2 n,o "bends" strong enough with n. This is definitely true on the expression 

= 3(1 + 2^) 2 " (8) 
derived in [6]. In the inertial range the dissipation contributions to (3) can be neglected. This 
does not mean that the even-order structure functions are not affected by the dissipation 
processes. The equation (3) is not closed and as a result the even -order moments are coupled 
to the dissipation contributions appearing in the equations for the odd-order moments. This 
will be discussed below. 

The equation (3) is a direct consequence of the Navier-Stokes equations. We will show in 
what follows that the ESS is consistent with (3). Let us, in accord with ESS (2), assume 
that S 2n = S 2n (S 2m ) where m is an arbitrary number. This assumption is non-trivial since, 
in principle, the moment S 2n can also depend on the displacement r and dissipation scale r\ 
(Reynolds number). Substituting this into (3) gives: 

dS 2n , _ (d - l)S 2n ,o ~ (d - l)(2n - 1)5^-2,2 + (2w - l)rV Xj2n -2 - (2w - l)vrD u ^ n -2 
dS 2m fi (d - l)S 2m ,o - (d - l)(2m - l)S 2m - 2t2 + (2m - l)rV x , 2m - 2 - (2m - l)vrD Uy2m -2 

(9) 

The relation (2) holds if the right side of (9) is equal to J^'°^° . Again, the relations (9) 
are exact everywhere as long as r/Lf — > 0. The mean field approximation, introduced in 
[6], is a statement that the pressure-gradient difference is expressible in term of a quadratic 
form of velocity differences. Since < Ap y (Au) 2 >=< Ap y (Av) 2 >= 0, we are left with: 

a(Au) 2 + b(Avf d / a \ 2 

Ap x = ^ '- ^ J - + c—(Au) 2 + ■ ■ ■ 10 

r ar 



The coefficients a, b and c etc are chosen so that Ap x = Ap x Au = Ap x Av = 0. We also 
have (see [6]): 

Ap y oc AuAv/r (11) 
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The equations (3) (9)-(ll) are not closed since we not have the relations coupling S 2ri: o with 
5*271-2, 2- We know that in the dissipation range, r/rj — > the functions S 2n fl oc 5 , 2n „ 2 ,2 and 
^2n,o = 2n, while in the inertial range the correlation functions are characterized by the non- 
trivial exponents (1). In principle, based on [6], we can easily write equations for S 2n -2,2- 
However, they involve the correlation functions S , 2 „_4 i 4 etc. 

Now we would like to ask the central question: consider a relatively low Reynolds number 
flow , so that the dissipation contributions to (3) cannot be neglected and the functions 
02n,o(O, ^) vary with the displacement r. What is the structure of the theory preserving (2) 
but strongly violating the inertial range scaling S 2nj0 oc r 52n -°? At the top of the dissipation 
range r/rj m 1 — 10 the scaling functions, violating the inertial range scaling are not small 
(see (1)). For 2m = 2 the equation (9) simplifies: 



dS 2nfi = (d - l)S 2nfl - (rf - l)(2n - l)S 2 n-2,2 + (2n - \)rV x ^ 2 - (2w - l)i/r£) U|2n _ 2 
dS 2fi (d-l)(S 2 , -S , 2 ) 

(12) 

where in incompressible, isotropic and homogeneous turbulence (d— 1) (So, 2 — 5 , 2j0 ) = — r ^f^ 
(see (3)). Both dissipation and pressure contributions do not appear in the denominator of 

(2n,0 

(12). The form of relation (12) tells us that the ESS S 2nfi = C 2ri:2 (S 2fi ) « 2 .o , exact in the 
inertial range if the relations (1) are valid, is possible only in an interval where the numerator 
in (12) is equal to -r^°. 

Substituting this into (12) and using the scaling form (1) gives: 



O0 2n fi _ dr _ <52n,0>J2n,0 T &n,o<p2n,o{x) dx 



dS2,0 ^ ^ 6,0^2,0 1 + 



(13) 



dr ' ' ' t,2nfl<t>2,o(x) dr 



i2n,0 
(2,0 



By assumption, S 2n = S 2n (S 2 ), subject to "boundary condition" S^o — C 2 n, 2 S 2 2 ' as 

(2n,0 

x — > oo. Solution to (13), satisfying these constraints, is: 2rij o oc (p 2 ^'° . Indeed, substituting 
this into the second equation (13), we are left with the differential equation, equivalent to 
the ESS (2) with the Reynolds-number-independent coefficient C 2rij2 . One can see that the 
ESS is the only universal solution to the equation (13), not involving any dependence on the 
Reynolds number (77). 
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Since (3) and (9) are a direct consequence of the equations of motion for velocity field, we 
conclude that the the ESS with non-trivial scaling exponents is consistent with the Navier- 
Stokes equations as long as the scaling assumption (1) is valid. 

The function (f>2,o( x ) can be readily self-consistently found from the well-known differential 
equation: 

S _ .8r + 6z/^ (14) 
dr 

The inertial range calculations [6] and both numerical and physical experiments [3] give 
£2,0 ~ 0.7 and <f>2,o(x) ~ 02,0 ~ 2.0 (Kolmogorov constant Ck ~ 1-6). Substituting the ESS 
expression S^q oc S 2:0 into (14) gives: 



6d(p 



2.0 



: =(0.8 - O.30j o )x 1 ^-660 2 ,o/x (15) 

where by definition of the dissipation scale vri~ 2+ & = 1 and £ — 1. Solution to this equation 
gives 02,0(^)5 g en tly approaching a 2) o ~ 2 as x — > 00. Noticible deviations from this constant 
value start at x « 30 — 50 (at x = 20 the function 2 ,o ~ 1-65 — 1.7). This equation was 
derived by Benzi et. al. (Ref.[2]). 

In the inertial range, where the dissipation contributions to (12) are negligible and where 
the ESS is sinonimous to the power laws, the numerator of (12) is equal to — r dS ^'° and the 
mean field approximation as exact as the power laws themselves. It is interesting that even 
if ^2n,o 7^ S 2n -2,2, the power laws solutions of (3) are still possible: in principle the mixed 
moments S 2n -2,2 can be cancelled by the corresponding pressure-gradients contributions to 
(10). 

To conclude: It follows directly from the Navier-Stokes equations that if the inertial range 
scaling exists, then: 

(d - l)S 2nfi ~{d- l)(2n - l)S 2n _ 2i2 + (2n - l)rV x , 2n - 2 = -r—^- (16) 

dr 

This expression means that the inertial range pressure contribution to this equation must 
be 0(S 2n fi) or 0(S 2n - 2y2 ). This proves the mean-field approximation (10). 
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It can be shown that in the interval x > 1, the direct dissipation contribution to (12) 
is small. Then, the mean -field approximation justifies the assumption S 2m fl = ^m.ot^o)- 
This leads to the ESS. 

However, the general statement, not related to a particular dynamical system, can be 
made: 1. if the scaling relation (1) with the 0(1) coefficients c n>m is valid; 2. if S n>0 = 

xi n,Q 

S n ,o{S m ,o) is independent on r and rj, then S nt0 = C ntm S^Q° . This relation means that there 
exist only one dominating (dynamically relevant) scaling function faj and all others can 
calculated in a simple way. 

Now we can discuss the cases where the ESS is violated. In a strongly sheared wall 
flow one can introduce two Reynolds numbers. The first one is Re = UL/u where L is the 
width of the channel (boundary layer, etc) and U is a characteristic (mean) velocity. The 
second one (Re = u^Ljv ) is based on the friction velocity ul = —v^\ W au- The dissipation 
rate £ = 0(^Re A jRe) is a weak function of the Reynolds number (dissipation ) scale. 
Thus, all structure functions, even if they can be written in a form (1), must involve the 
Re-dependent proportionality coefficients. This violates the assumptions leading to the ESS 
(2). Far enough from the wall, where £ ~ If 3 / L with the Re- independent proportionality 
coefficient one can expect the ESS to be valid. 

In some sheared flows 1/L S = j^/u = 0(1), the scaling functions also depend on y/L s fa 
1. In these regions the y/L s cannot be neglected and the simple derivation of the ESS (13) 
breaks down. The best example, illustraiting this point, is the Kolmogorov flow driven by 
the forcing function f = (0,cosL s x). There we expect the ESS to hold in the vicinities of 
zeros of the local strain rate d x U y oc sin(L s x) and break down near local maxima (minima) 
of the strain rate. These conclusions agree with the experimental findings [9], [10]. 
I am grateful to S.Kurien, K. R. Sreenivasan , A. Polyakov and M. Vergassola for their help 
and comments . 
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